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This contribution to the proceedings of the Cracow meeting on 'Applications of Random Matrix 
Theory' summarizes a series of studies, some old and others more recent on financial applications 
of Random Matrix Theory (RMT). We first review some early results in that field, with particu- 
lar emphasis on the applications of correlation cleaning to portfolio optimisation, and discuss the 
extension of the Marcenko-Pastur (MP) distribution to a non trivial 'true' underlying correlation 
matrix. We then present new results concerning different problems that arise in a financial con- 
text: (a) the generalisation of the MP result to the case of an empirical correlation matrix (ECM) 
constructed using exponential moving averages, for which we give a new elegant derivation (b) the 
specific dynamics of the 'market' eigenvalue and its associated eigenvector, which defines an interest- 
ing Ornstein-Uhlenbeck process on the unit sphere and (c) the problem of the dependence of ECM's 
on the observation frequency of the returns and its interpretation in terms of lagged cross-influences. 

I. PORTFOLIO THEORY: BASIC RESULTS 

Suppose one builds a portfolio of N assets with weight w, on the ith asset. The (daily) variance of the portfolio 
return is given by: 

R 2 = ^2 "'■ t7 - ( '■j G .i"'r (1) 

ij 

where <r| is the (daily) variance of asset i and Cy- is the correlation matrix. If one has predicted gains gt , then the 
expected gain of the portfolio is G = wigi. 

In order to measure and optimize the risk of this portfolio, one therefore has to come up with a reliable estimate of 
the correlation matrix Cy. This is difficult in general [lj,|2j since one has to determine of the order of N 2 /2 coefficients 
out of N time series of length T, and in general T is not much larger than N (for example, 4 years of data give 1000 
daily returns, and the typical size of a portfolio is several hundred stocks.) We will denote, in the following, q — N/T; 
an accurate determination of the true correlation matrix will require q <C 1. If r \ is the daily return of stock i at time 
t, the empirical variance of each stock is: 

^f>j) 2 , (2) 

and can be assumed for simplicity to be perfectly known (its relative mean square-error is (2 + k)/T, where k is the 
kurtosis of the stock, known to be typically k s» 3). In the above definition, we have, as usual, neglected the daily 
mean return, small compared to daily fluctuations. The empirical correlation matrix is obtained as: 



V rp 



1 T 

= 5>Hi x\ = r\j(Ji. (3) 



If T < N, E has rank T < N, and has N — T zero eigenvalues. Assume there is a "true" correlation matrix C 
from which past and future x\ are drawn. The risk of a portfolio constructed independently of the past realized x\ is 
faithfully measured by: 

( R b) = Wl<Jl ( x t x t) w j a j ~ X! "•-' 7 < ( "< ; rr ;"V (4) 

ijt ij 

Because the portfolio is not constructed using E, this estimate is unbiased and the relative mean square-error is small 
(~ 1/T). Otherwise, the w's would depend on the observed x's and, as we show now, the result can be very different. 

Problems indeed arise when one wants to estimate the risk of an optimized portfolio, resulting from a Markowitz 
optimization scheme, which gives the portfolio with maximum expected return for a given risk or equivalently, the 
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minimum risk for a given return (G): 

luij 9il ViCij 9j/<7j 

From now on, we drop a\ (which can always be absorbed in gi and wi). In matrix notation, one has: 



gT C -lg 



The question is to estimate the risk of this optimized portfolio, and in particular to understand the biases of different 
possible estimates. We define the following three quantities: 



• The "In-sample" risk, corresponding to the risk of the optimal portfolio over the period used to construct it 

_G 2 

gT E -l g 



R in = w^Ewb = „ T ^_, (7) 



• The "true" minimal risk, which is the risk of the optimized portfolio in the ideal world where C would be 
perfectly known: 



G 

#true = w C Cw C = j r _, (8) 
g J C x g 



• The "Out-of-sample" risk which is the risk of the portfolio constructed using E, but observed on the next period 
of time: 

2 rp ,g T E- 1 CE" 1 g 

#out = %Cw £ = G ( g T E -lg)2 (9) 

From the remark above, the result is expected to be the same (on average) computed with C or with E', the ECM 
corresponding to the second time period. Since E is a noisy estimator of C such that (E) = C, one can use a convexity 
argument for the inverse of positive definite matrices to show that in general: 

(g T E- 1 g) > g T C- 1 g (10) 

Hence for large matrices, for which the result is self-averaging: 

Rfn < -Rtrue- (H) 

By optimality, one clearly has that: 

fltnie - R ouf ( 12 ) 

These results show that the out-of-sample risk of an optimized portfolio is larger (and in practice, much larger, see 
Fig 1) than the in-sample risk, which itself is an underestimate of the true minimal risk. This is a general situation: 
using past returns to optimize a strategy always leads to over-optimistic results because the optimization adapts to 
the particular realization of the noise, and is unstable in time. In the case where the true correlation matrix is C = 1, 
one can show that Q: 

Q2 

#Le = -T- and R fn = ^rucV 1 ~<Z = ^outl 1 ~ S) (13) 

g g 

Only in the limit q — > will these quantities coincide, which is expected since in this case the measurement noise 
disappears. In the other limit q — > 1, the in-sample risk becomes zero since it becomes possible to find eigenvectors 
(portfolios) with zero eigenvalues (risk), simply due to the lack of data. 
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Risk 

FIG. 1: In sample (left curves) and out of sample (right curves) portfolio risk along the optimal curve in the risk return 
plane. Red lines are with the raw empirical matrix, blue lines with the cleaned matrix using RMT, showing how the risk 
underestimation can be reduced by matrix cleaning. From Q. 

II. MATRIX CLEANING AND RMT 

How should one 'clean' the empirical correlation matrix to avoid, as much as possible, such biases in the estimation 
of future risk? In order to get some general intuition on this problem, let us rewrite the Markowitz solution in terms 
of the eigenvalues Xk and eigenvectors V k of the correlation matrix: 

wnxY, V^»;;</ ; = * + £( A fc 1 - wtvfgj (w) 

kj kj 

The first term corresponds to the naive solution: one should invest proportionally to the expected gain (in units where 
<ji = 1). The correction term means that the weights of eigenvectors with A > 1 are suppressed, whereas the weights 
of eigenvectors with A < 1 are enhanced. Potentially, the optimal Markowitz solution allocates a very large weight to 
small eigenvalues, which may be entirely dominated by measurement noise and hence unstable. A very naive way to 
avoid this is to go for the naive solution w% oc <?j, but with the k* largest eigenvectors projected out: 

rvixgi- v iV*9i, (15) 

k<k*;j 

so that the portfolio is market neutral (the largest eigenvector correspond to a collective market mode, » 1/y/N) 
and sector neutral (other large eigenvectors contain sector moves). Since most of the volatility is contained in the 
market and sector modes, the above portfolio is already quite good in terms of risk. More elaborated ways aim at 
retaining all eigenvalues and eigenvectors, but after having somehow cleaned them. A well known cleaning corresponds 
to the so-called "shrinkage estimator" : the empirical matrix is shifted 'closer' to the identity matrix. This is a Bayesian 
method that assumes that the prior empirical matrix is the identity, again only justified if the market mode has been 
subtracted out. More explicitely: 

E c = orE + (1 - a)l so \ k c = 1 + a(\ k - 1), (16) 

where the subscript c corresponds to the cleaned objects. This method involves the parameter a which is undetermined, 
but somehow related to the expected signal to noise ratio. If the signal is large, a — > 1, and a — > if the noise is large. 
Another possible interpretation is through a constraint on the effective number of assets in the portfolio, defined as 
(^^luf) -1 Q. Constraining this number to be large (ie imposing a minimal level of diversification) is equivalent to 
choosing a small. 

Another route is eigenvalue cleaning, first suggested in 0, Q| , where one replaces all low lying eigenvalues with a 
unique value, and keeps all the high eigenvalues corresponding to meaningful economical information (sectors): 

\ k c = l-6 if k>k*, A*=A| if k < k*, (17) 

where k* is the number of meaningful number of sectors and 8 is chosen such that the trace of the correlation matrix 
is exactly preserved. How should k* be chosen? The idea developed in 0, Q is to use Random Matrix Theory to 
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determine the theoretical edge of the 'random' part of the eigenvalue distribution, and to fix k* such that A|, is close 
to this edge. 

What is then the spectrum of a random correlation matrix? The answer is known in several cases, due to the work 
of Marcenko and Pastur [f| and others We briefly recall the results and some elegant methods to derive them, 

with special emphasis on the problem of the largest eigenvalue, which we expand on below. Consider an empirical 
correlation matrix E of N assets using T data points, both very large, with q = N/T finite. Suppose that the true 
correlations are given by {x\x 3 t ,) = CijS tt i. This defines the Wishart ensemble 9]. In order to find the eigenvalue 
density, one introduces the resolvent: 



from which one gets: 



G(z) = ±Tr[(zI-E)- 1 ] (18) 



p(A) = lim-3(G(A-ie)). (19) 

c— >0 TT 



The simplest case is when C = I. Then, E is a sum of rotationally invariant matrices SEfj = {x\x\)jT. The trick 
in that case is to consider the so-called Blue function, inverse of G, i.e. such that B(G(z)) = z. The quantity 
R(x) = B{x) — 1/x is the 'R-transform' of G, and is known to be additive Since <5E* has one eigenvalue 

equal to q and N — 1 zero eigenvalues, one has: 



Inverting SGt(z) to first order in 1/N, 



and finally 



6B ^ = \ + W^) - B ^ = \ + ^y (21) 



r (r s (z + q-1)- V(z + g -l) 2 -4i^ 

Ge{z) = — , (22) 



which finally reproduces the well known Marcenko-Pastur (MP) result: 



pW _ ^- 2 ^-'> 2 (23) 

The case of a general C cannot be directly written as a sum of "Blue" functions, but one can get a solution using 
the Replica trick or by summing planar diagrams, which gives the following relation between resolvents: 0, 0i0 

zG E (z) = ZGc(Z) where Z= \ (24) 

1 + q{z&E{z) - 1) 

from which one can easily obtain p(X) numerically for any choice of C y|. [In fact, this result can even be obtained 
from the original Marcenko-Pastur paper by permuting the role of the appropriate matrices]. The above result does 
however not apply when C has isolated eigenvalues, and only describes continuous parts of the spectrum. For example, 
if one considers a matrix C with one large eigenvalue that is separated from the 'Marcenko-Pastur sea', the statistics 
of this isolated eigenvalue has recently been shown to be Gaussian p| ( see also below), with a width ~ T -1 / 2 , much 
smaller than the uncertainty on the bulk eigenvalues (~ q 1 ^ 2 )- A naive application of Eq. Ij24(l . on the other hand, 
would give birth to a 'mini-MP' distribution around the top eigenvalue. This would be the exact result only if the 
top eigenvalue of C had a degeneracy proportional to N. 

From the point of view of matrix cleaning, however, these results show that: (i) the expected edge of the bulk, 
that determines k* , obviously depends on the prior one has for C. The simplest case where C = I was investigated 
in particular in 0,^1, with the results shown in Fig 2. Other choices are however possible and could lead to an 
improved cleaning algorithm; (ii) the uncertainty on large eigenvalues is much less than that on the bulk eigenvalues, 
meaning that the bulk needs a lot of shrinkage, but the bigger eigenvalues less so - at variance with the naive shrinkage 
procedure explained above. An alternative route may consist in using the 'power mapping' method proposed by Guhr 
|l4j or clustering methods [Tlj . 
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FIG. 2: Empirical eigenvalue density for 406 stocks from the S&P 500, and fit using the MP distribution. Note (i) the presence 
of one very large eigenvalue, corresponding to the market mode (see section IV) and (ii) the MP fit reveals systematic deviations, 
suggesting a non trivial structure of the true correlation matrix, even after sector modes have been accounted for (see |3.ll6j|1. 



III. EWMA EMPIRICAL CORRELATION MATRICES 



Consider now the case where C = I but where the Empirical matrix is computed using an exponentially weighted 
moving average (EWMA). More precisely: 



t-i 



E tj =e £ {X-if-t'xtxi (25) 



V — — oo 



with e = 1/T. Such an estimate is standard practice in finance. Now, as an ensemble 25y satisfies Eij = (1 — e)£V 



exjxj. We again invert the resolvent of JEt to find the elementary R-function 



SB t (x) = -+R t (x) with R t ( x ) = —1 (26) 

x N{1 — qx) 

Using the scaling properties of G(z) we find for R(x): 

R aC (x) = aRc(ax). (27) 

This allows one to write: 

R E (x) = R (l - t)E {x) + R t (x) = (1 - q/N)R E ((l - q/N)x) + (28) 

To first order in 1/N, one now gets: 

R(x) + xR'ix) + — g— = — > R(x) = - l0g(1 ~ qx) . (29) 
1 — qx qx 

Going back to the resolvent to find the density, we finally get the result first obtained in [l7|: 

p(X) = -3G(A) where G(X) solves XqG = q — logfl - qG) (30) 

This solution is compared to the standard MP distribution in Fig 3. 

Another nice properties of Blue functions is that they can be used to find the edges of the eigenvalue spectrum 
(A±). One has: 12] 

X±=B(x±) where B'(x±) = (31) 
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FIG. 3: Spectrum of the exponentially weighted random matrix with q = Ne = 1/2 and the spectrum of the standard Wishart 
random matrix with q = N/T — 1/3.45, chosen to have the same upper edge. From [t^I . 



In the case at hand, by evaluating B(x) when B'(x) — we can write directly an equation whose solutions are the 
spectrum edges (A±) 

A± =log(A ± )+ 9 +l (32) 

When q is zero, the spectrum is a S in 1 as expected. But as the noise increases (or the characteristic time decreases) 
the lower edge approach zero very quickly as A_ ~ exp(— q). Although there are no exact zero eigenvalues for EWMA 
matrices, the smallest eigenvalue is very close to zero. 



IV. DYNAMICS OF THE TOP EIGENVALUE AND EIGENVECTOR 



As mentioned above, it is well known that financial covariance matrices are such that the largest eigenvalue is well 
separated from the 'bulk', where all other eigenvalues reside. The financial interpretation of this large eigenvalue is 
the so-called 'market mode': in a first approximation, all stocks move together, up or down. One can state this more 
precisely in the context of the one factor model, where the zth stock return at time t is written as: 

= + 4 (33) 

where the market mode <pt is common to all stocks through their market exposure (3 l and the e\ are idiosyncratic 
noises, uncorrelated from stock to stock. Within such a model, the covariance matrix reads: 

d^frpjol + olSii. (34) 

When all <7j's are equal, this matrix is easily diagonalized; for N stocks, its largest eigenvalue is A = • Pj) a j> + (j2 
and is of order N, and all the other N — 1 eigenvalues A Q are equal to a 2 . The largest eigenvalue corresponds to the 
eigenvector More generally, the largest eigenvalue Ao, normalized by the average square volatility of the stocks, 
can be seen as a proxy for the average interstock correlation. 

A natural question, of great importance for portfolio management, or dispersion trading (option strategies based 
on the implied average correlation between stocks), is whether A and the j3's are stable in time. Of course, the 
largest eigenvalue and eigenvector of the empirical correlation matrix will be, as discussed at length above, affected 
by measurement noise. Can one make predictions about the fluctuations of both the largest eigenvalue and the 
corresponding eigenvector induced by measurement noise? This would help separating a true evolution in time of the 
average stock correlation and of the market exposure of each stock from one simply related to measurement noise. 
We shall see that such a decomposition seems indeed possible in the limit where Ao 3> A a . 

We will assume, as in the previous section, that the covariance matrix is measured through an exponential moving 
average of the returns. This means that the matrix E evolves in time as: 



E ijt t = {l-e)E ijt t- 1 + er\r{. 



(35) 
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The true covariance matrix Cy = (rV J ) is assumed to be time independent - this will give us our benchmark 
hypothesis - with its largest eigenvalue A associated to the normalized eigenvector ^> . In this section we deal with 
covariance matrices instead of correlation matrices for simplicity, but most results should carry over to this latter case 
as well. 

Denoting as Aot the largest eigenvalue of E t associated to ipot, standard perturbation theory, valid for e <C 1, gives: 

Aot = (1 - e)A t-i + e(V>ot-i|C|V>ot-i) + e (^ot-i\VtHot-i) , (36) 
with r\ij = rV J — (rV J ). We will suppose for simplicity that the returns are Gaussian, yielding: 

(VtjVke) = CikCjg + CuCjk- (37) 
In the limit where Ao becomes much larger than all other eigenvectors, the above equation simplifies to: 

Aot « (1 - e) Aot-i + e cos 2 t _iA o [1 + &] , (38) 

where cos9 t = (V'otl^o) and £t is a random noise term of mean zero and variance equal to 2. This noise becomes 
Gaussian in the limit of large matrices, leading to a Langevin equation for Ao: 

^ = e(cos 2 6A -A ) + ecos 2 6>6- (39) 
at 

We have neglected in the above equation a deterministic term equal to e sin 2 9Ai , which will turn out to be a factor 
(Ai/Ao) 2 smaller than the terms retained in Eq. (|38J) . 

We now need to work out an equation for the projection of the instantaneous eigenvector ipot on the true eigenvector 
^o- This can again be done using perturbation theory, which gives, in braket notation: 



Wot) = Wot-i) + e) — \ip. 



,„ Aot_l — Aat-l 
a^tO 



at-1 



_ , , ■ , r t r t \ipot-i) (^ot-i|rtrt|^ot-i) | . , . 

~ \Wot-i)+e e l^ot-i), (40) 

Aot-l -Aot-i 

where we have used the fact that the basis of eigenvectors is complete. It is clear by inspection that the correction 

term is orthogonal to |-0ot l) s so that \ipot) is still, to first order in e, normalized. Let us now decompose the matrix 

r t r t into its average part C and the fluctuations r], and first focus on the former contribution. Projecting the above 
equation on < leads to the deterministic part of the evolution equation for cosf?t: 

cos 6>t » cos#t_i + ecos# f _i- — ecos 3 0t-\- — — , (41) 

Aot-i Aot-i 

where we have neglected the contribution of the small eigenvalues compared to Ao, which is again a factor (Ai/Ao) 2 
smaller. In the continuous time limit e — > 0, this equation can be rewritten as: 

and describes a convergence of the angle 9 towards or 7r - clearly, 'i'o and — ^o are equivalent. It is the noise term 
•q that will randomly push the instantaneous eigenvector away from its ideal position, and gives to 9 a non-trivial 
probability distribution. Our task is therefore to compute the statistics of the noise, which again becomes Gaussian 
for large matrices, so that we only need to compute its variance. Writing \ipot) — cos9t\^o) + sin^tl^i*)) where \^u) 
is in the degenerate eigenspace corresponding to small eigenvalues Ai , and using Eq. I|37|) , we find that the noise term 
£t to be added to Eq. (|4*2"|l is given by: 

e 2 

(Ct 2 ) » tt- [2A^ cos 2 t sin 2 9 t + A Ax cos 2 29 t ] , (43) 

where we have kept the second term because it becomes the dominant source of noise when 0^0, but neglected 
a term in A 2 . The eigenvector ipo therefore undergoes an Ornstein-Uhlenbeck like motion on the unit sphere. One 
should also note that the two sources of noise £t and are not independent. Rather, one has, neglecting A 2 terms: 

2(6Ct) « A cos 2 t sin20 t - Ai sin49 t (44) 
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FIG. 4: Variogram of the top eigenvalue of the correlation matrix Ao (top curve) and (log)-variogram of the top eigenvector, 
In {{^ot+rlV'ot}}, as a function of the time lag r, for US stocks, and for an exponential moving average parameter e = 1/30 days. 
We show for comparison (bottom curve) the Ornstein-Uhlenbeck behaviour expected for both quantities in the case of a time 
independent underlying correlation structure. 

In the continuous time limit, we therefore have two coupled Langevin equations for the top eigenvalue Xot and the 
deflection angle t . In the limit Ai <C A , the stationary solution for the angle can be computed to be: 



P{6)=N 



l+cos20(l-fc) 

Ai' 
Ao 



-i l/4e 



1 -cos 20(1 



(45) 



As expected, this distribution is invariant when 9 — > n — 9, since — V&o is also a top eigenvector. In the limit Ai Ao, 
one sees that the distribution becomes peaked around 9 = and it. For small 9, the distribution becomes Gaussian: 



P{9) 



1 



: exp 



2e^ 



(46) 



leading to (cos 2 9) w 1 — eAi/2A The angle 9 is less and less fluctuating as e — > (as expected) but also as Ai/A — > 0: 
a large separation of eigenvalues leads to a well determined top eigenvector. In this limit, the distribution of Ao also 
becomes Gaussian (as expected from general results |l3j ) and one finds, to leading order: 



(Ao) w A - eAi/2; ((SX ) 2 ) 



(47) 



Therefore, we have shown that in the limit of large averaging time and one large top eigenvalue (a situation approx- 
imately realized for financial markets), the deviation from the true top eigenvalue <5Ao and the deviation angle 9 are 
independent Gaussian variables (the correlation between them indeed becomes zero as can be seen using Eq. (I44|l in 
that limit, both following Ornstein-Uhlenbeck processes. 

From these results one directly obtains the variogram of the top eigenvalue as: 



([A f+T -A t ] 2 ) = 2e (1 - cxp(-er)) . 



(48) 



This is the result expected in the absence of a 'true' dynamical evolution of the structure of the matrix. From Fig. 
4, one sees that there is clearly an additional contribution, hinting at a real evolution of the strength of the average 
correlation with time. One can also work out the average overlap of the top eigenvector with itself as a function of 
time lag, E((tpo t+T \ipo t )). Writing again \if>ot) = cos0t|^o) + sin0{|^it)) one sees that an equation for the evolution of 
the transverse component \^u) is a priori also needed. It is easy to show, following the same method as above, that 
the evolution of the angle <pt made by the component with a fixed direction is a free diffusion with a noise term of 
order eAi/A . Therefore, on the time scale e _1 over which 9 t evolves, \^u) can be considered to be quasi-constant, 
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leading to: 

((V'ot+rlV'ot)) » £(cos(# t - e t+T )) » 1 - e-i(l - exp(-er)). (49) 

Ao 

Any significant deviation from the above law would, again, indicate a true evolution of the market structure. Again, 
Fig. 4, provides some evidence of such an evolution, although weaker than that of the top eigenvalue Aq. 



V. FREQUENCY DEPENDENT CORRELATION MATRICES 

The very definition of the correlation matrix a priori depends on the time scale over which returns are measured. 
The return on time r for stock i is defined as: ri tT (t) = \npi(t + r) — hxpi(t), where pi(t) is the price at time t. The 
correlation matrix is then defined as: 

Q , (r) = M>Vif (50) 

A relatively well known effect is that the average inter-stock correlation grows with the observation time scale - this 
is the so-called Epps effect [l8|, [19J . For example, for a collection of stocks from the FTSE index, one finds, in the 
period 1994-2003: 

(Qyy(5'))«0.06; w0.19; (Qyy(ld)) « 0.29 (51) 

Besides the change of the average correlation level, there is also a change of structure of the correlation matrix: the 
full eigenvalue distribution distribution changes with r. A trivial effect is that by increasing the observation frequency 
one also increases the number of observations; the parameter q defined above decreases and the noise band is expected 
to shrink. This, at first sight, appears to be a nice way to get rid of the observation noise in the correlation matrix 
(see [20] for a related discussion). Unfortunately, the problem (or the interesting effect, depending on the standpoint) 
is that this is accompanied by a true modification of the correlations, for which we will provide a model below. In 
particular one observes the emergence of a larger number of eigenvalues leaking out from the bulk of the eigenvalue 
spectrum (and corresponding to 'sectors') as the time scale r increases. This effect was also noted by Mantegna pi): 
the structure of the minimal spanning tree constructed from the correlation matrix evolves from a 'star like' structure 
for small t's (several minutes), meaning that only the market mode is present, to a fully diversified tree at the scale 
of a day. Pictorially, the market appears as an embryo which progressively forms and differentiates with time. 

The aim of this section is to introduce a simple model of lagged cross-influences that allows one to rationalize the 
mechanism leading to such an evolution of the correlation matrix. Suppose that the return of stock i at time t is 
influenced in a causal way by return of stock j at all previous times t' < t. The most general linear model for this 
reads: 

n,i{t) = &(*) + 2 / dt'Kij(t - t') rj>1 (t') (&(*)&(«')) = DiSySQ - (52) 
3 J -°° 

Here r = 1 is the shortest time scale - say a few seconds. The kernel is in general non-symmetric and describes 
how the return of stock j affects, on average, that of stock i a certain time later. We will define the lagged correlation 
Cij{t-t')h T . 

Ci^t-t 1 ) = (ri,i(t)r jtl (t')}. (53) 

This matrix is, for t 7^ t', not symmetric; however, one has obviously Cij(t — t') — Cji(t' — t). These lagged correlations 
were already studied in |22|. Going to Fourier space, one finds the Fourier transform of the covariance matrix 
Cij(uj) =Cji{-w)\ 

Cy(w) = £(1 - K{u:))7 k \l - K(-u>))j£D k (54) 
k 

where K(ui) is the Fourier transform of K(t) with by convention K(r < 0) = 0. When cross-correlations are small, 
which is justified provided the Ti t x(t) corresponds to residual returns, where the market has been subtracted, the 
relation between Cij and Kij becomes quite simple and reads, for r > 0: 

Cy(r) = DjKijir). (55) 
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FIG. 5: Typical self-influence kernel Ka(E) in Laplace space and fit with the Laplace transforms of the sum of two exponentials. 
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FIG. 6: Typical cross-influence kernels Kij(E) for three pairs of stocks, and fit with the Laplace transforms of the sum of two 
exponentials. Note that the influence amplitudes have different signs, even for the same pair of stock, depending on the time 
scale. 



This equation allows one, in principle, to determine Kij(r) from the empirical observation of the lagged correlation 
matrix. Suppose for simplicity that the influence kernel takes the form Kij(r) — K^e~ TijT , then Kij(oj) — K®j/(iu> + 
^ij). In this model, the primary object is the influence matrix K which has a much richer structure than the correlation 
matrix: each element defines an influence strength K° and an synchronisation time F . In fact, as shown in Figs. 
5 and 6, fitting the empirical data requires that Kij is parameterized by a sum of at least two exponentials, one 
corresponding to a time scale of minutes, and a second one of hours; sometimes the influence strength corresponding 
to these two time scales have opposite signs. Pooling together the parameters corresponding to different pairs of 
stocks, we find, as might have been expected, that strongly coupled stocks (large K°) have short synchronisation 
times r . 

Coming back to the observation that the correlation matrix is frequency dependent, one should note that the scale 
dependent correlation matrix Cij(r) is related to Cij(u>) by: 

Cy( T ) = { r i,Trj, T )c = / dto S 2 (u>T)Ci : j(uj) (56) 
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where S(.) is the form factor (i.e. Fourier transform of the window used to define returns on scale r, for example a 
flat window in the simplest case). Therefore, for r small one finds that residuals are uncorrelated (i.e. the correlation 
matrix has no structure beyond the market mode): 

C ij (r->Q)KiD i 6 ij , (57) 

whereas on long time scales the full correlation develops as: 

p oo 

Cy(T->oo)»A*0+ / dTiDjKiji^+DiKjiir)]. (58) 
Jo 

The emergence of correlation structure therefore reveals the lagged cross-influences in the market. Note that on long 
time scales, small K 0, s can be counterbalanced by large synchronisation times and lead to significant correlations 
between 'weakly coupled' stocks. 

We believe that a more systematic empirical study of the influence matrix Kij and the way it should be cleaned, 
in the spirit of the discussion in section II, is worth investigating in details. 

We want to thank Pedro Fonseca and Boris Schlittgen for many discussions on the issues addressed in sections IV and V, 
and Szilard Pafka and Imre Kondor for sharing the results on the EWMA matrices given in section III. We also thank Gerard 
Ben Arous and Jack Silverstein for several clarifying discussions. We also thank the organisers of the meeting in Cracow for 
inviting us and for making the conference a success. 
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